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SUMMARY 


The  indicial  lift  and.  pitching-moment  coefficients  are  derived,  for 
flat-plate  triangular  wings  traveling  at  supersonic  speeds.  The  coeffi¬ 
cients  are  determined  for  angle— of -attack  distributions  corresponding  to 
sinking  wings  and  to  pitching  wings.  The  wing  with  supersonic  edges  is 
completely  analyzed,  and  the  wing  with  subsonic  edges  is  partially 
analyzed,  the  solution  in  this  case  being  completed  for  very  narrow 
wings  by  an  application  of  slender  wing  theory.  In  the  case  of  the 
supersonic  edges,  a  comparison  is  made  with  known  two-dimensional 
results  and  also  with  the  results  for  the  same  triangular  wing  in 
reversed  flow. 


INTRODUCTION 


The  wing  of  triangular  plan  form  has  received  considerable  atten¬ 
tion  in  the  steady-state  theory  of  three-dimensional  wings  in  a  super¬ 
sonic  stream.  The  purpose  of  the  present  report  is  to  determine  the 
aerodynamic  characteristics  of  a  triangular  wing  in  supersonic  unsteady 
motion. 

There  are  several  simple  types  of  unsteady  motion  on  which  the 
analysis  can  be  based.  The  so-called  indicial  motion,  in  which  the 
velocity  undergoes  a  discontinuous  change  at  t*=0,  will  be  considered 
here.  (See  also  references  1  and  2.)  It  is  possible  to  conceive  the 
physical  situation  in  two  slightly  different  ways.  For  one,  it  can  be 
supposed  that  the  wing  has  been  traveling  at  the  constant  velocity  VD 
for  an  infinitely  long  time  and  then,  at  t‘=0,  starts  suddenly  to  sink 
without  pitching  motion  (or  to  pitch  without  sinking)  while  maintaining 
the  forward  velocity  V0.  On  the  other  hand,  the  wing  may  be  considered 
to  be  at  rest  in  still  air  until  at  t*=0  it  starts  suddenly  either  to 
sink  or  to  pitch  and,  at  the  same  instant,  attains  the  forward  velocity 
V0.  The  latter  physical  picture  will  be  used  in  this  report.  Problems 
of  unsteady  motion  can  also  be  approached  with  the  initial  assumption 
that  the  velocity  potential  depends  harmonically  on  the  time.  (See 
reference  3.)  These  two  approaches  are  quite  compatible  in  that  they 
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can  be  related  through  the  use  of  superposition  methods  (Duhamel *  s 
integral ,  Fourier fs  integral)  of  the  operational  calculus. 

The  question  of  whether  the  coordinate  system  should  move  with  the 
wing  or  remain  fixed  is  also  of  some  importance.  (See  reference  4  for 
a  discussion.)  The  latter  alternative,  that  is,  where  the  wing  moves 
away  from  the  coordinate  system,  was  chosen  for  this  report  because  the 
velocity  potential  cp  in  this  case  satisfies  the  wave  equation 

^  xx  +  ^yy  +  ^  zz  “  ^"2  ^tft 1  =  ® 

where  x,y, z  are  Cartesian  coordinates,  t*  is  time,  and  aQ  is  the 
speed  of  sound  in  the  undisturbed  medium.  The  fact  that  the  equation 
has  this  form  is  helpful  in  establishing  analogs  between  steady  and 
nonsteady  motions,  and  these  analogs  are  of  considerable  help  in  the 
solution  of  certain  problems. 

The  boundary  conditions  to  be  considered  correspond  to  the  problem 
of  the  flat  sinking  wing  (angle— of— attack  distribution  uniform  over  the 
plan  form)  and  to  the  flat  pitching  wing  (angle  of  attack  varies 
linearly  with  chordwise  distance)  in  indicial  motion.  First,  the  load¬ 
ing  on  a  flat  triangular  wing  with  supersonic  edges  undergoing  an 
indicial  sinking  motion  is  determined.  Then  a  simplified  method  is 
developed  whereby  total  lift  and  pitching-moment  coefficient  for  the 
wing  with  supersonic  edges  may  be  obtained.  These  quantities  are 
determined  as  functions  of  time,  for  both  sinking  and  pitching  wings. 

Lastly,  the  triangular  wing  with  subsonic  edges  is  partially 
analyzed,  and  an  approximate  method  for  very  slender  wings  is  used  to 
complete  the  determination  of  loading.  The  analog  method,  mentioned 
previously,  is  here  of  great  value.  Lift  and  pitching-moment  coeffi¬ 
cients  for  the  sinking  and  pitching  slender  triangular  wing  are 
determined. 


LIST  OF  IMPORTANT  SYMBOLS 


speed  of  sound  in  the  free  stream 


cQ  root  chord  of  triangular  wing 

/  lift 

Ct  lift  coefficient  l  - — 

L  ^§Oo  V02S 


) 


indicial  lift  coefficient  due  to  angle— of-attack  change  (without 
pitching)  (cl.  *  ^U) 
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£p 

% 

a 


s 

S 
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indicial  lift  coefficient  due  to  pitching  on  a  wing  rotating  about 
its  leading  edge  or  apex 


pitching-moment  coefficient,  positive  when  trailing  edge  tends  to 
sink:  relative  to  leading  edge 

indicial  pitching- moment  coefficient  due  to  angle— of-attack  change 
(without  pitching)  measured  about  the  leading  edge  or  apex 


indicial  pitching-moment  coefficient  due  to  pitching  measured  about 

^m1 


the  leading  edge  or  apex 


ciV 


cotangent  of  sweep  angle  (cot  A) 
free— stream  Mach  number  fe) 


°  V 


o  5=0 


Jfe  I**  ay 

2as  J_ g  q.o 

mMoVo  f  d 
2  8  s2  4  9o 


loading  coefficient  (pressure  on  the  lower  surface  minus  pressure 
on  the  upper  surface  divided  by  free-stream  dynamic  pressure) 


free— stream  dynamic  pressure 


dimensionless  rate  of  pitching 


local  semi  span  of  wing 

wing  area 

time 
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*0 

U,V,W 


Vo 


*>7»  z 


^n.p. 


t 

Co 

perturbation  velocity  components  in  the  x,y,z  directions, 
respectively 

free— stream  velocity 

Cartesian  coordinates 

distance  of  center  of  pressure  from  wing  apex 


a  angle  of  attack  (angle  between  flight  path  and  plane  of  wing), 

radians 

0  yii-Mo2i 


Pe 


1 


m 


Mq  4  n/i+eP" 


0 

m 

Q 


A 


wing  angle  of  pitch,  relative  to  initial  attitude,  positive  when 
trailing  edge  lies  below  leading  edge 


wing  rate  of  pitch,  positive  when  trailing  edge  is  sinking 


relative  to  leading  edge 


angle  of  sweep  of  leading  edge ,  positive  for  sweepback 


pQ  free-stream  density 


T 


t 

S 


chord  lengths  traveled 


or  Mot0 


<P 

£q> 


perturbation  velocity  potential 

jump  in  potential  across  the  z=0  plane 
[<p(x,y,0+)  -  cp(x,y,0  )] 


A 

A 


« 


Subscripts 

component  taken  normal  to  the  leading  edge 

positive  side  of  the  z=0  plane,  or  upper  surface  of  a  wing 
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THE  PROBLEM  AND  THE  NATURE 
OF  ITS  SOLUTION 


Consider  a  wing  situated  in  still  air  and  a  Cartesian  coordinate 
system  with  origin  on  the  leading  edge  at  the  point  of  wing  symmetry. 

At  time  equal  to  zero,  the  wing  starts  impulsively  to  move  in  a  straight 
line  with  constant  velocity  away  from  the  coordinate  system  which 
remains  fixed  relative  to  the  still  air  at  infinity.  The  load  distri¬ 
bution  on  a  wing  undergoing  such  a  motion  is  called  the  indicial  loading. 
Similarly,  the  forces  and  moments  which  are  based  on  this  loading  are 
given  the  adjective  indicial.  The  partial  differential  equation  that 
is  satisfied  by  the  velocity  potential  T  for  such  a  motion  can  be 
written 

^tt  “  ^xx  ^yy  -<Pzz  =  0  (***) 

which  is  the  normalized  form  of  the  wave  equation.  In  equation  (l), 
x,y,  and  z  are  distances:  x  measured  chordwise,  y  spanwise,  and  z 
vertically,  and  t  is  equal  to  a0t!  where  aQ  is  the  free— stream 
ft  speed  of  sound  and  t*  is  time. 

The  boundary  conditions  to  which  equation  (l)  is  subject  are 
^  dependent  on  the  wing  shape  and  motion.  Adopting  the  assumptions  of 

thin-airfoil  theory,  which  are  consistent  with  the  assumptions  already 
used  in  obtaining  equation  (l),  it  can  be  assumed  that  the  slope  of  the 
wing  surface,  in  the  direction  of  motion,  at  any  place  and  time  is 
given  by  the  ratio  of  the  vertical  velocity  component  in  the  z=0  plane 
to  the  wing 1 s  forward  velocity  component.  For  a  flat  plate,  then,  the 
following  conditions  are  to  be  satisfied:- 

1.  The  vertical  velocity  cp^  is  a  linear  function  of  x, 
the  coefficients  of  which  depend  upon  the  angle  of  attack  and 
rate  of  pitch,  over  the  portion  of  the  xy  plane  occupied  by 
the  wing  at  any  given  time. 

2.  No  perturbations  exist  at  infinity. 

3*  There  are  no  discontinuities  in  the  velocity  potential 
except  over  the  region  occupied  by  the  wing  and  its  vortex  wake. 

The  problem  has  now  been  expressed  as  one  of  finding,  for  pre¬ 
scribed  boundary  conditions,  a  solution  to  the  wave  equation.  It  is 
*  often  desirable  to  express  the  solution  in  terms  of  the  loading  coeffi¬ 

cient  rather  than  the  potential  function  or  velocity  components.  This 
coefficient  can  be  written  in  its  linearized  form  as 


6 


NACA  TN  2387 


Ap  2  dACp  /2 

<lo  VpMp  -5F 

The  solution  of  the  problem  without  further  restrictions  has,  in 
general,  not  been  obtained,  although  it  can  be  shown  that  the  boundary 
conditions  for  any  wing  in  unsteady  motion  can  be  satisfied  by  a 
suitable  superposition  of  sources  and  doublets  (reference  4).  However, 
the  solution  for  wings  with  all  supersonic  edges  can  be  written  as  a 
double  integral  of  sources  having  intensities  determined  by  the  local 
slope  of  the  wing.  Hence,  for  a  flat  surface  at  constant  angle  of 
attack  and  not  pitching. 


cp(x,y)  =  -§“  f  f  F 


(3) 


where  SQ  is  the  outline  of  the  region  of  sources  which  can, 

a, 


given  time,  affect  the  point  x,y 
determined,  r  equals  J (x— Xi )2 
attack  of  the  wing.  The 
form  and  a  discussion  of 


at  a 

at  which  the  potential  is  being 
is  the  angle  of 


+  (y-Ji)  ;  and  a 
area  Sa  has  been  termed  the  acoustic  plan 
its  significance  is  given  in  reference  4. 


The  solution  for  the  triangular  wing  with  subsonic  edges  can  be 
obtained  in  certain  regions,  but  in  others  the  problem  reverts  to  the 
solution  of  a  double  integral  equation  involving  time  and  the  two  sur¬ 
face  dimensions  of  the  wing.  If  the  triangular  wing  is  slender,  an 
approximate  method  for  finding  the  pressure  over  the  entire  wing  can  be 
used.  This  method  is  to  neglect  the  streamwise  velocity  gradients  in 
comparison  with  the  gradients  in  the  plane  normal  to  the  free  stream 
and  also  in  comparison  with  the  term  cp-tt*  There  results  for  the 
partial  differential  equation  governing  the  flow  field  (equation  (l)) 
the  wave  equation  of  one  lower  dimension,  namely, 

-  v  -  - 0  (l* 


As  will  be  developed  later,  the  boundary  conditions  become  the  same  as 
for  a  rectangular  flat  plate  of  very  low  aspect  ratio  inclined  at  angle 
of  attack  to  a  free  stream  with  a  Mach  number  equal  to  \/27  This 
analogy  with  the  steady-state  lifting— surface  problem  is  useful  since 
solutions  to  the  latter  problem  have  been  obtained. 

The  indicial  lift  and  pitching-moment  coefficients  will  be  deter¬ 
mined  for  wings  with  two  different  vertical  velocity  distributions. 

The  first  is  the  case  in  which  the  boundary  condition  is 


A 

* 
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where  a  is  the  angle  of  attack,  and  wu  is  constant  over  the  plan 
form.  This  case  is  referred  to  as  that  of  the  sinking  wing.  The  other 
"boundary— value  problem  considered  is  the  one  corresponding  to  a  flat 
wing  pitching  at  a  constant  rate  about  its  leading  edge  or  apex.  The 
boundary  conditions  have  the  form 

WU  =  ^2^=0  =  (x+MOt) 

where  9  is  the  constant  rate  of  pitch  and  (x  +  Mot)  is  the  distance 
aft  the  leading  edge  or  apex  of  the  wing.  The  pitching  velocity  Q 
is  considered  positive  when  the  trailing  edge  of  the  wing  sinks  relative 
to  the  axis  of  rotation. 


The  difference  between  a 
and  0  is  illustrated  in  the 
sketch,  where  for  clearness  the 
motion  is  shown  to  be  oscillatory 
rather  than  indie ial.  The  angle 
of  attack  a  is  the  angle 
between  the  flat  wing  surface  and 
the  line  tangent  to  the  flight 
path  of  the  leading  edge  or  apex 
of  the  wing.  The  angle  0  is 
the  angle  between  the  flat  wing 
surface  and  the  horizontal  (see 
part  (a)  of  the  sketch).  Part  (b) 
of  the  sketch  shows  a  wing  under¬ 
going  a  sinusoidal  angle-of-attack 
variation  with  a  zero  angle  of 
pitch  throughout  the  motion. 

Part  (c)  shows  a  wing  undergoing 
a  sinusoidal  angle-of-pitch  varia¬ 
tion  taken  about  the  le'ading  edge, 
with  the  angle  of  attack  remain¬ 
ing  zero. 

Although  the  lift  and 
pitching-moment  coefficients  are 
given  only  for  the  two  types  of 
motion  separately,  it  is  possible, 
according  to  the  linear  theory, 
to  combine  them  to  simulate  an 


Direction  of  wing  motion 

- > 


0=9=0 


arbitrary  indicial  maneuver  consisting  of  both  sinking  anri  pitching 
motions.  For  such  a  case,  the  boundary  condition  becomes 


wu  =  — VQa  -9  (x  +  ^t) 
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If,  furthermore,  the  desired  maneuver  is  not  indicial,  the  lift  and 
pitching-moment  coefficients  may  he  determined  by  use  of  Duhamel’s 
integral: 


CL  =  -jfr  f  (t,-T4)  a(T’ )  +  ^'(t’-r*  )q.(  t*)]  dT* 

cm  =  jT  [^’(t'-  t*  )a(  r*  )  +  ^’(t*-  t*  )  q(  t*  )]  dr' 

vhere  a(t’)  and  q(t*)  are  the  arbitrary  motions  ^q  =  ,  C^’, 

Cp^1,  and  are  the  indicial  aerodynamic  coefficients,  and  the  primes 

on  the  coefficients  indicate  that  the  pitching  motion  is  about,  and  the 
pitching  moments  are  referred  to,  the  leading  edge  or  apex. 


The  position  of  the  axis  of  pitching  motion,  and  of  the  axis  to 
which  pitching-moment  coefficient  is  referred,  is  of  importance.  In 
this  report,  these  axes  coincide  in  an  axis  normal  to  the  root  chord  of 
the  wing  and  passing  through  the  wing  leading  edge  or  apex.  However, 
it  is  often  desired  to  transfer  the  pitching  motion  and  moment  calcula¬ 
tion  to  other  axes,  and  the  formulas  for  such  a  transformation  will  be 
given  here.  Let  the  pitching  motion  refer  to  an  axis  lying  a  distance 
acQ  back  of  the  leading  edge  or  apex,  and  let  the  pitching-moment 
coefficients  refer  to  an  axis  bc0  aft  the  leading  edge  or  apex.  The 
necessary  transformation  formulas  are 


Cmalb  "  C™a'  +  b  CLa 

CLn  I  =  CT  '  -  a  Ct 
1j<1  Lq  -nx 

a 

« 

Cmq  Iq  b  =  C®4*  +  b  °Lq’  “  a  “  ab  cLa 

where  a  prime  on  a  quantity  indicates  that  the  pitching  motion  is  about, 
and  the  pitching  moments  are  measured  about,  an  axis  through  the  leading 
edge  or  apex  of  the  wing.  #The  quantity  q  is  the  dimensionless  rate 
of  pitching,  equal  to  cQ  9/VQ.  The  subscripts  a  and  b  mean, 
respectively,  that  the  quantity  in  question  refers  to  a  pitching  motion 
about  an  axis  at  a  distance  acc  aft  the  leading  edge  or  apex,  or  that 
the  pitching  moments  are  measured  about  an  axis  bc0  aft  the  leading 
edge  or  apex. 


NACA  TN  2387 


9 


WING-  WITH  SUPERSONIC  EDGES 
INDICIAL  LOADING  FOR  SINKING  WING 
Analysis 


The  seven  regions,—  The  analytic  expression  for  the  indicial  load¬ 
ing  over  the  triangular  wing  has  a  different  form  in  each  of  seven 
regions.  These  regions  are  determined  *by  the  positions  of  the  various 
wave  fronts  relative  to  the  wing  plan  form  (fig.  1).  For  t<0  the 
wing  is  motionless,  its  leading  edge  lying  along  lines  represented  hy 
the  dashed  lines  in  figure  1.  At  t=0  the  wing  starts  suddenly  to 
move.,  and  for  t>0,  travels  forward  at  a  constant  speed  VQ.  After  a 
certain  time  t  has  elapsed.,  the  wing  has  traveled  to  a  new  position, 
also  shown  in  the  figure.  In  this  same  interval  of  time,  pressure 
impulses  have  traveled  out  in  spherical  waves  from  every  point  of  the 
region  which  the  wing  has  occupied.  The  trace  on  the  wing  of  the 
sphere  starting  from  the  wing  apex  at  t=0  forms  the  external  boundary 
of  region  7.  The  area  outside  this  circle  and  within  the  traces  of  the 
cylindrical  waves  (the  envelopes  of  the  spherical  waves)  generated  hy 
the  leading  edges  at  t=0  forms  region  4.  Region  5  is  formed  hy  the 
overlapping  of  these  cylindrical  waves,  and  the  solution  for  loading 
within  it  can  he  found  hy  a  suitable  superposition  of  the  solutions  for 
regions  3  and  4.  Region  1  lies  between  the  cylinder  trace  on  the  wing 
and  the  leading-edge  position  at  time  t;  the  loading  in  this  region 
cannot  he  affected  hy  the  manner  in  which  the  wing  started  its  motion 
since  it  lies  outside  the  starting  cylindrical  waves.  .  Hence,  the 
loading  in  region  1  is  the  same  as  that  on  a  swept  wing  flying  at  a 
steady  supersonic  speed.  The  solution  in  region  2  can  also  he  obtained 
from  steady— state  lifting— surface  theory,  hut,  whereas  In  region  1  the 
field  is  two-dimensional  (i.e.,  invariant  with  distance  measured  parallel 
to  the  leading  edge),  in  region  2  the  field  is  conical.  Region  6  is 
formed  hy  the  overlapping  of  regions  2  and  4.  Finally,  region  3  is  that 
area  completely  unaffected  hy  waves  from  the  wing  edges.  In  the  follow¬ 
ing  subdivisions  the  analysis  of  each  of  the  separate  regions  will  he 
discussed. 

Region  1:  The  loading  in  region  1  of  figure  1  is  equal  to  the 
loading  on  a  two-dimensional  flat  plate  moving  at  a  constant  velocity 
given  hy  the  component  of  stream  velocity  normal  to  the  leading  edge 
of  the  triangular  wing.  Since  this  component  is  supersonic,  the 
loading  is  of  the  Ackeret  type  and  is  given  hy 

(  ^ 

V  9n/x  pn 

But  since  Vn  =  VQ  cos  A  where  A  is  the  angle  of  sweep  (see  fig.  2), 
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=  q.0  cos2a 

Cjj  =  a  sec  A 
^  =  Mq  cos  A 

Pn  =A2  cos2A-1 


and 


4a 


sec2A 


Finally,  if  ctn  A  =  m 


( 


Ap\ 

<lo  / 


4xm 


1  Jfl  -  1 


where  P  =a/ Mq2  -  1 . 


(5) 


Region  2:  The  steady-state  loading  on  a  triangular  wing  with 
supersonic  edges  has  been  given  by  several  authors  (see,  for  conven¬ 
ience,  reference  5)  so  the  expression  for  the  loading  in  region  2  can 
be  written  immediately  for  the  coordinate  system  shown  in  figure  2  as 

/  Ap  \  _  4am  -  {  tt  +  arc  sin  fsMs&Sj - 

VW2  7t/?5U>-  «m(I+Mot)-,] 


arc  sin 


pgmy+(x+Mnt )  \ 
p[M(x+Mot  )+y]  J 


(6) 


Region  3:  Since  region  3  is  unaffected  by  the  edges  of  the  wing 
the  solution  for  the  loading  therein  can  be  written  as  in  reference  6 


(7) 


Region  4:  The  solution  for  loading  in  region  4  can  be  obtained 
from  consideration  of  a  two-dimensional  wing  starting  from  rest  and 
moving  with  velocity  Vn  normal  to  its  leading  edge.  This  problem 
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has  teen  treated  in  reference  6  and  the  solution  written  there  can  he 
written  for  the  right-hand  side  of  figure  1  as 


1! 

r 

— j 

MnXn+t  Jj  Mn2— 1 

arc  cos  -  +  — 77 - 

f  —  +  arc  sin  —  ^ 

K  % X 

TT|3n  L 

xn+^t  Mn 

v 2  V  J 

where  the  notation,  as  defined  by  the  sketch,  is 


xn  =  x  cos  A  —  y  sin  A 
yn  =  x  sin  A  +  y  cos  A 


and  since 

ctn  A  =  m,  sin  A  =  —  ^  -  , 

*/  1+m2 
m 

cos  A  =  —==. 

*/  1+m2 


Then 


mx-y 
*/  1+m2 


x+my 
•J  1+m2 


The  equation  for  loading  now  becomes,  in  the  coordinate  system  of 
figure  2, 


4am 


it*/  p2]# 


-1  - 


arc  cos 


mMo(mx  —  |y|  )  +  tCl+m2) 
Vl+m2  (mx  —  |y  |  +  mMot ) 


*/  p2m2-l  f Tt  mx  —  |  y|  ^ 

— - -  1  A  +  arc  sin  — -  1 

^o  \2  tj  1+m2  / 


(8) 


Region  The  solution  for  loading  in  region  5  can  he  obtained 
by  superposition  of  the  solutions  for  regions  3  and  If  the  solu¬ 
tions  for  the  two  sides  of  region  4  (obtained  from  equation  8)  are 
added,  the  result  gives  a  value  of  wu  of  twice  the  required  amount 
in  region  5*  as  well  as  undesirable  pressures  off  the  wing.  However, 
subtraction  from  this  sum  of  the  solution  for  region  3  (equation  7) 
reduces  the  downwash  wu  to  the  proper  value,  and  also  cancels  the 
excess  pressures.  The  resulting  expression  can  be  written 


12 


KACA  TN  2387 


4am 

«  j  P2™2— 1 


mM^mx-y)  +  (l+m?)t 
arc  cos  — — — — -  + 

v'l+m2  (mx— y  +  mMot) 


arc  cos 


inMo  (mx+y)  +  (l+n? )  t 
a/  1+m?  (mx  +  y  +  mMot) 


(  arc  sin  -JSg  -  +  arc  sin  .  ^ 
“Mq  V  y'J+P'  t  /TtP  t 


(9) 


Region  6:  The  loading  in  region  6  can  also  be  calculated  by 
superposition.  To  find  the  loading  in  this  case,  add  the  solutions 
for  regions  2  and  4  (equations  (6)  and  (8))  and  subtract  the  solution 
for  region  1  (equation  (5))*  There  results 


4am 

« *J  02m2_l 


■j^arc  sin 


Pgmy  -  (x+Mpt) 
PL’m(x+Mot)  -y] 


—  arc  sin 


Pgmy  +  (x+Mpt ) f 
ptm^x+Mpt  )+y] 


mMo(mx  —  |  y|  )  +  t(l+mg) 
(mx  —  |  y  |  +  mMot  )V  1+m2 


+ 


^  h  f  —  +  arc  sin 
eiMq  \  2 


nnc  -  lyl  ^  X 

Wl+nf2  '  * 


(10) 


Region  7?  The  solution  for  the 
obtained  by  means  of  equation  (3). 


loading  in  region  7  can  be 

The  analysis  used  in  finding  the 
solution  in  this  region  is 
not  difficult  but  the  algebra 
is  rather  involved.  It  is 
useful  at  this  point  to 
introduce  polar  coordinates 
(see  sketch)  such  that 

x— xi  =  r  cos  9 

y-y  1  =  r  sin  9 
dxjdyi  =  r  dr  d0  (11) 

From  equation  (ll),  equa¬ 
tion  (3)  can  be  written  in 
the  form 

cp  =  f  f  drdR  (12) 

2it  js^J 
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The  acoustic  plan  form  for  points  in  region  7  is  the  region 
hounded  by  three  curves  as  indicated  in  the  sketch.  The  arc  between 
0i  and  02  is  determined  by  eliminating  T  between  the  equations 

r2  =  (t  -  T  f 

(the  equation  for  the  inverse  sound  waves)  and  the  equation  for  the 
left  leading  edge 

7i  =  -m(xi  +  MoT) 

The  arc  between  %  and  0i  is  found  by  determining  the  acoustic 
intersection  (eliminating  T)  of  the  right  leading  edge  with  the 
inverse  sound  wave;  and,  finally,  the  arc  between  9 2  and  £3  is 
given  by  the  equation  r=t .  The  equations  of  these  arcs  can  be 
written  in  polar  coordinates  as 
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- - siaSiite  ;8lS95te 

m  cos  0+mMo+sin  0 


r  =  t;  9z  <  9  <  93 


m(x+Mot )-y  . 

m  cos  0+idMq— sin  9  ’ 


0  <  9  <  9X 
9q<9  <2rt 


Using  these  expressions,  equation  (12)  reduces  to 

cp  =  1°^  r92__m(x+M0t)+j  d0  + 

2n  Jg  m  cos  0+mMo+sin  9 


td0  + 


IsE  PdX  m(x+Mpt )— y  dQ 

2n  Jg  m  cos  9  +mMo— sin  9 

3 

and  taking  the  partial  derivative  with  respect  to  t  (to  determine 
the  loading  according  to  equation  (2))  one  finds1 


^  j  _  2a  m  j  _ d0 _  +  _2a  P 

<lo4  n  mMq  +  m  cos  0+sin  9  jiMq 


“  ^ 9  mMq  +  m  cos  9— sin  9 


In  evaluating  this  equation,  the  following  integral  is  used: 


for  — Jt  <  0  <  jt 


mMo  +  m  cos  0  ±  sin  0  »J  p2m2-l 


=  arc  tan  — .(Mg  ' 


l^m2-! 


1The  limits  02,  02,  and  03  are  all  functions  of  t  hut  in  moving 

the  partial  derivative  through  the  integral  sign  the  terms  involving 
vOi/dt,  and  all  cancel  one  another* 
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Since  equation  (l4)  is  valid  only  in  the  interval  -it  <.  0  ^  it,  care 
must  te  exercised  in  applying  it  "because  the  angle  93  may  be  greater 
than  rr  (as  in  the  preceding  sketch).  In  case  jt  <  9 3,  it  is  con¬ 
venient  to  introduce  the  angle  03f  =  03-2rt.  The  expression  for 
Ap/q0  can  ^hen  be  written  in  two  forms,  according  as  0Q  is  less 
than  or  greater  than 

for  y  >  0,  03  <n 

it  +  arc  tan  H.  _ 

J  P2n?-1 


In 


4am 


/q2™2 


* 


arc  tan  tan 


m  (Mq-1  )  tan  ( fl3  /2  )-l  ] 
arc  tan  - .  - 

J  62m2-l  J 


+  ,  x 

*Mo  v  3  2 ; 


(15a) 


for  y  >  0,  )t  <  03 
( 

V'IqA  it  Vp2®2-!  L 


arc  tan 


m(MQ-l)tan  (gi/2)-l 
V  p2]]!2— 1 


m(M0-l)tan  (0i/2)+l  ,  _  +_  m(Mo-l)tan  (92/2)+l 

arc  tan  -  ..  ■  -  — - -  +  arc  tan  - == - 

J  P2!!!2-!  yp2™2-! 


arc  tan 


m(Mo— l)tan  (03’/2)— 1 


+  5?1  (03’“®2  +  2it)  (15b) 

*M0  3 


where 


0i  =  arc  cos 

(0  <  0x  <  it) 


-Mo  y2  +  (x+Mot) 

y2  +  (x+MqI f 


02  =  arc  cos 


m(y+mx)  -  J~{ l+m?)t2  -  (y+mx)2 


(l+n?)t 


(0  <  02  <  it) 


1 6 
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@3  =  arc  cos 


-m(y-mx) 


—  J  (l+mF)t2  —  (y-mx)2 


(l+m2)  t 


03»=  03  -  2jt  (if  It  <  03) 


* 


The  limitation  on  03  can  be  given  both  an  analytic  and  geometric 
interpretation.  Thus,  equation  (15(a))  applies  for  0  <  m(x+t)<  y 
and  equation  (15(b))  applies  for  0  <  y  <  m(x+t).  These  regions  are 
shown  in  the  accompanying  sketch.  Because  of  the  geometrical  symmetry 
about  the  x  axis,  equations  (15)  suffice  for  the  determination  of 
loading  throughout  region  7» 


Discussion  of  Results 


g 


Plots  of  the  load  distribution  on  the  sinking  triangular  wing  with 
supersonic  edges  are  shown  in  figure  3 >  and  an  isometric  drawing  of  the 


NACA  TN  2387 


17 


loading  on  the  right  panel  appears 
in  the  sketch.  The  positions  of 
the  spanwise  sections  were  chosen 
so  that  each  of  the  regions  1 
through  7  is  represented.  It  is 
to  be  noted  that  the  results  for 
region  7  show  no  unusual  charac¬ 
teristics,  but  fit  in  well  with 
those  for  the  adjoining  regions. 

In  general,  the  distribution  is 
similar  to  the  steady— state  load¬ 
ing  on  a  triangular  wing. 


HTOICIAL  LIFT  AMD 
PITCHING  MOMENT 
Analysis 

Methods  of  solution.—  If  a 


detailed  knowledge  of  the  load 

distribution  is  not  required,  but  onl;y  the  total  values  of  lift  and 
moment  need  be  known,  analyses  much  simpler  than  the  one  presented  in 
the  previous  section  can  be  employed.  These  methods,  however,  not  only 
require  the  edges  of  the  wing  to  be  supersonic,  but  also  require  the 
trailing  edge  to  be  straight  and  normal  to  the  free— stream  direction. 
One  such  simple  method  has  been  presented  in  reference  4.  It  involves 
the  integration  over  the  three-dimensional  plan  form  of  longitudinal 
strips,  or  elements,  which  carry  the  two-dimensional  values  of  loading 
as  a  function  of  time.  The  results  presented  in  the  present  section 
can  be  derived  by  this  method  as  well  as  by  the  nfcthod  to  be  developed 
next. 


Consider  again  equation  (l)  and  integrate  each  term  with  respect 
to  y  between  the  limits  minus  and  plus  infinity.2  There  results  the 
equation 


If  7  =  y^(x,z,t)  and  y  =  yr(x,z,t)  are  the  equations  of  the  Mach 
waves  streaming  back  from  the  leading  edges  on  the  left  and  right  sides 
of  the  wing,  respectively,  then,  since  cp  is  continuous  across  these 

2  — —  —  —  —  ■ 

The  basic  idea  for  this  solution  was  given  by  Prof.  P.  A.  Lager strom 

in  his  lectures  at  the  California  Institute  of  Technology. 
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waves  but  cpx,  cpy,  cpz,  and  <Pt 


are  not. 


Syr  dyz 
dx  r  8x 


u2 


d2<P 

dx2 


<iy 


where  ur  and  uj  are  the  values  of  u  on  the  interior  faces  of  the 
right  and  left  Mach  waves,  respectively.,  and 


Values  of  the  terms  involving  w  and  are  similar  to  those  involv¬ 

ing  u  so  that  finally,  if 


then 


(16) 


2  ^ 

8  $ 

dx2  +  8z2 


f'^Jl 

\  8x 


uz 


+ 


Sy2 


K 


* 


4 


The  terms  enclosed  within  the  brackets  in  the  last  equation  combine 
so  that  each  bracketed  quantity  is  zero.  For  the  case  of  interest  here* 
this  is  not  difficult  to  show.  Consider,  for  example,  the  right  wedge. 
Then,  since  the  equation  of  the  wedge  is 


yr  = 


/  2  2 

— z  v  p  m  — 1  +  mx  +  ml^t 


and  the  value  of  the  potential  is  the  steady— state  two-dimensional  value 
given  by  the  expression 


— w, 


m(x+Mot)-y-z  J 
a/P2  v?-i 


X 
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3It  is  not  necessary  to  perform  a  direct  calculation  in  order  to  prove 
the  above  result  for  arbitrary  plan  forms*  The  bracketed  terms 
represent  the  directional  derivative  of  the  velocity  potential  taken 
along  the  so-called  "cononrial"  of  the  foremost  disturbance  surface. 
Since  cp  is  constant  on  the  surface,,  and  since  the  conormal  lies 
along  the  surface.,  the  bracketed  terms  are  zero. 
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for  a  flat  triangular  wing  with  a  constant  rate  of  pitch  0  about  its 
apex  (positive  9  produces  a  downward  motion  of  the  trailing  edge 
relative  to  the  leading  edge). 


These  boundary  conditions  are  exactly  like  those  studied  in  steady- 
state  supersonic  wing  theory.  In  fact,  the  lifting-surface  analog 
(shown  on  the  sketch)  is  a  wing  tip  of  specified  camber  in  a  supersonic 
free  stream  having  a  Mach  number  equal  to  The  solution  for  the 

potential  in  the  plane  of  the  wing  for  this  problem  can  therefore  be 
witten  immediately  as 


where  a  is  the  portion  of  the  area  on  the  shaded  surface  in  the  sketch 
lying  ahead  of  the  forecone  traces  given  by  the  equation  (t— ti)2  = 
(x-xi)2.  Using  equation  (2)  for  the  loading  coefficient,  and  introducing 
the  following  notation  for  the  average  spanwise  loading 

P°  =  2fia  f  ^dy  (sinkinS  wing)  (20a) 

Pi  =  ?  f  — dy  (pitching  wing)  (20h) 

2s*  9  u-s  <10 
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where 


s  =  mCx+Mot) 


it  is  found  that 

P  (  P  (n+Mpti)  dxidti 

0  ns  8t  Jg  J  yj t-ti)2  -  (x-xi)2 

Pi  ...  4m2  8  P  P  (xi+Moti)2  Ati 
*  1  ‘  ns2  St  Ja  J  /(t_tl)2  _  (x_Xl)2 


(21) 


(22) 


(23) 


Wing  vlth  constant  angle  of  attack  (sinking  wing).—  The  solution 

for  the  average  load  on  a  flat,  supersonic-edged,  triangular  wing 
starting  from  rest  at  t=0  and  flying  at  a  constant  speed  and  angle  of 
attack  is  given  hy  equation  (22).  With  the  transformations 


X-Xi  =  X2 
t— t  j_  =  tg 


this  becomes 


Po 


4m  a 

ns  St 


dtgdxg 


|  -  xg  -  Mptg 

V  tg2  -  Xg2 


This  integral  can  be  evaluated  and  gives  for  x<t  (region  A) 


P 


o 


=  4 


(24a) 


for  — t  <  x  <  t 


(region  B) 


a/^-x2  + 
x+MQt 


Mo 

arc  cos 


t+MoX 

X+Mgt 


and  for  x  <  — ’ t 


(region  C) 


(24b) 


p 


(24c) 


where  the  regions  are  shown  in  the  preceding  sketch.  Equations  (24)  can 
also  be  obtained  by  integrating  the  equations  for  the  loading  given  in 


22 


NACA  TN  2387 


the  preceding  section.  These  integrations  were  carried  out  (in  some 
regions  numerically)  and  the  results  were  found  to  agree  with  those  of 
the  present  analysis. 

It  is  now  possible  to  write  the  equation  for  the  indicial  lift  and 
pitching  moment  for  a  sinking  wing: 


2 

S  Mo 


pC0-M.Qt 

“Mot 


m(x+Mot)  P0  dx 


2 

S  Cq  Mo 


'C0-MOt' 


m(x+MQt)2  P0  dx 


(25) 

(26) 


where  S  is  the  wing  area  (equal  to  mc02)  and  the  prime  indicates  that 
the  pitching  moment  is  measured  about  the  apex,  the  positive  moment 
being  one  which  causes  the  trailing  edge  to  sink  relative  to  the  apex. 

Combining  equations  (24)  and  (25)  one  finds  for  the  first  interval 
shown  in  the  sketch 


a 


( 


Interval  no. 
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°L^  = 


a  Moco2  ^  >4 it  P 


,-t 


UMq 


(x+Mot)  dx  + 


(x+Mq^ ) 


-t 


x+Mot 


+  arc  cos 


(-f 


Mr>  t+MoX 

arc  cos  — -r^- 
3  x+MQt 


dx  + 


c0-4*ot 


4(x+Mot)  dx 


This  equation  integrates  to  give,  if  t0  =  t/c0, 

for  0  <  tD  <  — —  (first  interval) 

°  “  Mo+1 


CT  =-=-(l+lt2 

La  Mo 

Similarly  for  -----  tD  <  — ^ -  (second  interval) 


Mq+1 


Mq— 1 


Cj  =  -4—  (l  +  —  tQ2N)arc  cos  ^ + 

“  "Mo  \  2  0  J  tQ 


*3 


arc  cos  (M0— tQ32 )  + 


2  4t02  -  (1  -  Mot0); 


and  for  <  tQ  (third  interval) 

11q  ± 


Pl 


a  p 


In  the  same  manner  the  values  for  C^1  in  the  various  intervals 
determined  by  combining  equations  (24  j  and  (2 6).  There  results 


(27a) 


(27b) 


(27c) 
can  be 


2k 
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for  0  <  <  -  (first  interval) 

~  0  “  Mq+1 


V  - 


(28a) 


_ i _  <  <  — —  (second  interval) 

Mq+1  0  "  mq— 1 


C  *  - —  4  (84^t0-Jfe8t08-2t0a)  yt0z-(l^otor  + 

“a  «M0  Llo 


=r  (2+MQt03)  arc  cos  —  +  31,0  003  (Mq- 


(28b) 


and  for  — —  <  tQ  (third  interval) 


(28c) 


Numerical  results  for  a  Mach  number  of  2  will  be  presented  in  the  dis¬ 
cussion  section. 

Wing  with  linear  angle-of-attack  variation  (pitching  wing).—  The 

solution  for  the  average  load  on  a  flat,  supersonic— edged,  triangular 
wing  flying  at  a  constant  speed  and  pitching  at  a  uniform  rate  9 
about  its  apex  is  given  by  equation  (23  )•  With  the  transformation 

x— xi  =  xa 


t— t  x  =  t*2 


this  equation  becomes 


p  _  isf.  A 

1  ns2  8t 


-  -  x2  -  Mots 

75* 


dt2  dxg 
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and  the  evaluation  of  this  gives  (for  the  regions  defined  for  the 
sinking  wing) 

for  x  <  t  (region  A) 


p1  =  4  1  +  i.  (  — L_ 

L  2  \  x+Mot 


(29a) 


— t  <  x  <  t  (region  B) 


1  t2+(x+Mot)‘ 

2  (x+Mot  )2 


arc  cos 


i) + 


&  arc  003  *  1 

13  x+M0t  2  (x+Mot)2 


(29b) 


and  for  x  < -t  (region  C) 


Pi  = 


(29c) 


The  equations  for  lift  and  pitching  moment  on  a  pitching  wing  can 
"be  obtained  from  the  equations 


=  CT  .  -  _iL_ 
c09\  8  ScQ  Mq 

V„  J 


c0-M0t 


(x+Mot)2  Px  dx 


5[T  =  (V 

.  YrsJ 


Sc°  Mo  <4l0t 


Co-V 


(x+Mot)  P1  dx 


where  the  primes  indicate  the  wing  is  pitching  about  and  the  moments  are 
measured  about  the  leading  edge* 

.  +A  combination  of  equations  (29)  and  (30)  gives  for  the  lift  coeffi- 
c  lent 


for  0  <  t„  < 


o  —  ■  (first  interval) 
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v  ■  C1  + t to*  ■  M°t°3) 


(32a) 


•fnr  1  <  t  2  — —  (second,  interval) 

Mq+1  Mo-1 


%  jtMq 


(l  tc>2  “  +  |) 


Mo^Or-l 

arc  cos  — ; - * 


Mo 

3P 


arc  cos 


(Mq— p2t0)  +  ( |  to2  +  |  ^^o2-  M°t°  + 

jj)  ^t02-(i-Mot0r  ]  (32t) 


and  for  -lr  <  tQ  (third  interval) 

Mq-1 

V  -  8/3P 

Similarly  a  combination  of  equations  (29)  and  (31)  yields,  for  the 
pitching  moment  about  the  apex,  the  results 

for  0  <  tn  <  — i—  (first  interval) 

Mq+1 


(32c) 


P  I  —  _ 

H  “  Mo' 


1  +  to2  -  |  to4  (1+1Mo2 ) 


(33a) 


for  1  <  to  <  rrlr  (second  interval) 

Mq+1  0  Mq-1 


P  t  =  _  ±. ... 

tcMq 
2Mo 


2(l+t02)  -|t04  (1+4*0 


arc  cos 


Mot0-l 


— -2.  arc  cos  (Mq— p  tQ)  + 

itp 


jl_ 

12 


42-22MotQ+ (2Mq2+3 )t02+ (2Mq2 +13 )Moto 


^/tQ2— (l— Mot0 )  f  (33^) 
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and  for  - 

Mq-1 


<  tD  (third  interval) 


(33c) 


Numerical  results  for  a  Mach  number  of  2  are  presented  in  the 
next  section. 


Discussion  of  Results 


Figure  4  shows  the  variation  of  CL^,  C^',  and  the  location  of 
xc.p.  with  tq,  the  number  of  chord  lengths  traveled  (  tq  =  V0t * /cG  = 
Mc,t/c0  =  Mq^)  for  a  free-stream  Mach  number  Mq  equal  to  2.  The 
values  of  and  C^’  are  in  agreement  with  those  given  in  refer¬ 

ences  1  and  2.  For  the  purpose  of  comparison,  similar  curves  are  shown 
for  a  two-dimensional  wing  having  a  chord  cQ  equal  to  the  root  chord 
of  the  triangular  wing,  and  for  the  triangular  wing  in  reversed  flow. 
The  material  necessary  for  the  calculation  of  the  latter  curves  was 
presented  in  reference  6.  Several  interesting  conclusions  can  be  drawn 
from  these  results. 

First,  notice  that  the  total  indicial  lift  on  the  triangular  sink¬ 
ing  wing  is  the  same  at  every  instant  as  that  on  the  same  wing  in 
reversed  flow  (both  wings  of  course  having  started  with  the  same  velo¬ 
city  at  the  same  time),  and  that  the  value  of  this  lift  is  the  ««Tna  as 
the  total  indicial  lift  on  the  two-dimensional  wing  only  at  the  begin¬ 
ning  of  the  motion  and  again  when  the  steady  state  has  been  attained 
(fig.  4(a)).  Such  a  result  for  reversed  flow  is  not  true  for  the 
pitching  moment  (fig.  4(b))  and  center  of  pressure,  and  it  can  be  shown 
that  it  is  true  for  the  total  lift  only  when  the  wing  is  a  flat  lifting 
surface  with  supersonic  edges.4 

Second,  notice  that,  since  all  of  the  characteristics  for  the 
triangular  wings  are  independent  of  the  angle  of  sweep,  they  are  valid 
for  any  unyawed  triangular  wing  flying  at  a  Mach  number  equal  to  2 
and  having  supersonic  leading  edges. 

Third,  it  is  apparent  from  figure  4(a)  that  the  transition  of  the 
total  indicial  lift  from  its  initial  to  its  final  value  is  less  abrupt 
than  that  transition  for  the  two-dimensional  wing,  nnd  finally  the 
movement  of  the  center  of  pressure  on  the  sinking  wing  of  all  types  is 
seen  from  figure  4(c)  to  be  small. 


R.  T.  Jones  has  shown,  in  an  unpublished  work,  that  the  buildup  of 
total  indicial  drag  on  symmetrical  nonlifting  wings  is  the  same  for 
all  types  of  plan  forms  in  forward  and  reversed  flow. 
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Finally ,  it  can  "be  seen  from  equations  (27a)  and  (27c)  that 
the  initial  and  final  values  of  Ct  depend  on  1/Mq  and  1/P,  respec- 

CL 

tively.  As  the  Mach  number  Mg  is  increased,  therefore,  the  variation 
dies  out  since  p  and  Mq  become  nearly  equal.  The  same  renark 
applies  to  all  the  other  coefficients. 

Figure  5  presents  the  values  of  C]^  ' ,  C^^*,  an^-  ^-c .p. /co 

wings  pitching  about  the  foremost  extremity  of  their  plan  form  (i.e., 
leading  edge  or  apex).  Again  the  results  are  presented  in  terms  of 
chord  lengths  traveled  tq  for  a  free-stream  Mach  number  equal  to  2. 

For  these  wings  it  is  apparent  that  the  reverse-flow  theorem  does 
not  apply  even  to  total  lift.  The  results  for  the  triangular  wing  are 
still  independent  of  the  angle  of  sweep,  however,  and  the  movement  of 
the  center  of  pressure  is  again  slight. 


WING  WITH  SUBSONIC  EDGES 
INDICIA!  LOADING  FOR  SINKING  WING 
Analysis 


The  six  regions.-  As  in  the  study  of  supersonic-edged  triangular 
wings,  there  are  also  in  the  case  of  triangular  wings  with  subsonic 
leading  edges  various  regions  in  which  the  analytical  form  of  the  load¬ 
ing  equation  is  different.  Figure  6  shows  the  regions  into  which  the 
subsonic-edged  triangular  wing  can  be  most  conveniently  divided.  Most 
of  these  regions  have  counterparts  on  the  supersonic-edged  wing  shown 
in  figure  1. 

To  begin  with,  region  6  lies  within  the  spherical  wave  which 
started  at  t=0  from  the  wing  apex.  Region  1  is  within  the  cylindrical 
wave  which  was  started  at  t=0  by  one  edge  of  the  wing,  but  outside 
the  wave  started  by  the  other  edge •  Region  4  is  the  area  formed  by  the 
overlapping  of  the  two  cylindrical  waves  from  the  opposite  edges,  but 
outside  the  region  influenced  by  the  reflection  of  one  of  these  waves 
on  the  opposite  edge  (secondary  wave  fronts  shown  in  the  figure). 

Region  5  is  the  area  between  regions  4  and  6  where  the  flow  is  influ¬ 
enced  by  secondary  (and  higher— order )  wave  reflections.  Finally, 
regions  2  and  3  are  similar  to  regions  2  and  3  in  the  supersonic-edged 
case;  region  2  being  that  uninfluenced  by  the  starting  phenomena  and 
therefore  having  a  loading  already  at  Its  steady— state  value,  and 
region  3  being  that  which  is  unaffected  by  the  disturbances  emanating 
from  the  edges. 
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Region  1:  The  solution  for  the  load  distribution  in  region  1 

is  the  same  as  that  for  a  t wo-d ime n 3 i onal  wing  starting  suddenly  from 
rest  and  moving  with  a  steady  subsonic  velocity  vn.  A  solution  to 
the  latter  problem  for  the  initial  part  of  the  motion  is  presented  in 
reference  7«  In  terms  of  the  normal  components  of  velocity  and  dis¬ 
tance^  therefore,  the  loading  coefficient  for  the  right-hand  side  of 
figure  6  can  be  written  immediately: 


8v0 

*VnMn 


(  Jk_  / ~ 

\1+Mn  V  Mnt+xn 


+  arctM/SHi) 

V  t-xn  / 


The  equations  which  relate  the  normal  components  to  those  in  the 
free— stream  direction  have  already  been  given  in  the  section  on 
region  4  of  the  supersonic-edged  wing.  Use  of  these  relations  leads 
to  the  following  expression  for  loading  in  region  1  (in  the  coordinate 
system  of  fig.  7) 


Ap 


8a 


mMo 


t  y  1+m2  +  [  y  1  — 


mx 


qQ>l  tfMo  \mMo+  J±+tx£  V  m  M0t  -  |y|  + 


mx 


arc  tan 


mMpt  —  [y  |  +  mx 
t  V l+m2  +  |y|  —  mx 


(3*0 


Region  2:  The  loading  on  region  2,  being  the  steady— state 

loading  on  a  triangular  wing  with  subsonic  edges,  is  well  known.  The 
solution  for  region  2  of  figure  7  is  therefore  given  by  (see,  for 
convenience,  reference  5) 

( 


ap 

Q. 


O' 2 


4am2  (x+M0t ) 
m?  (x+Mot  )2— y2  E 


(35) 


where  E  is  the  complete  elliptic  integral  of  the  second  kind  with 
modulus  *J  1— @2m2 


Region  3:  The  loading  in  region  3  follows  from  reference  6  and 
is 


f  ApA  4a 
\  "  % 


(36) 
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Region  4:  The  loading  in  region  4  of  figure  7  is  calculated  by 
superposition.  Just  as  the  solution  for  region  5  of  the  wing  with 
supersonic  edges  was  obtained.  The  solution  in  region  4  is  the  sum 
of  the  solutions  for  the  right  and  left  halves  of  region  1,  minus  the 
result  for  region  3»  Thus 

( 


/mM^+y+mx 

—  - -  + 

t  1+mF— ; y-mx 


;v  l+n^+y-mx 
mMot— y+mx 


Ap^  -  ^  (  ft  Vi+mP-3Mmc 

^0/4  ^  _  mMo+  aJT+vP  \V  mMot+y+mx 


Regions  5  and  6:  In  these 
regions  the  exact  solution  for 
the  loading  has  not  "been  deter¬ 
mined.  As  was  shown  in  refer¬ 
ence  such  solutions  would 
require  the  solution  of  a  three- 
dimensional  elliptic— type  partial 
differential  equation.  In  this 
report  a  later  section  will  con¬ 
tain  an  approximate  solution  for 
these  regions. 

Discussion  of  Results 


An  isometric  drawing  of  the 
load  distribution^  for  the 
regions  in  which  it  is  known,,  is 
shown  in  the  sketch.  Comparing 
the  results  for  the  loading  on 
this  wing  to  the  one  with  super¬ 
sonic  edges  (fig.  3)*  tt  is 


apparent  that  the  principal  difference  in  the  two  distributions  is  in 
the  behavior  around  the  leading  edges;  the  loading  being  finite  at  the 


supersonic  edge,,  whereas  it  becomes  infinite  at  the  subsonic  edge.  In 
view  of  the  known  steady— state  results  this  difference  was  to  be 


expected.  Elsewhere  the  loadings  are  quite  similar. 
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The  results  presented  in  equations  (34)  through  (37)  will  next  he 
examined  in  a  different  light.  Choose  a  given  spanwise  section  on  the 
wing  and  watch  this  section  as  time  progresses  from  t=0.  This  amounts 
to  fixing  the  axis  on  the  body  and  can  he  accomplished  simply  by  using 
the  quantity  s  introduced  in  equation  (21), 

s  =  m(x  +  Mot) 

It  is  clear  that  s  is  the  semispan  of  a  given  spanwise  section,  and 
that  if  equations  (3*0  through  (37)  are  written  in  terms  of  s,  y,  and 
t,  for  a  fixed  s  they  represent  the  variation  of  loading  on  a  given 
section  as  time  progresses. 

If  the  notation  is  further  simplified  by  introducing  the  parameter 
Pe  where 


P  =  -  ■ 

r q  . . .  , 

mMo+  V 1+m2 


(38) 


equations  (34)  through  (37)  can  be  written  in  the  following  way: 


8a 

*Mo 


(aM.  JiK  \  B  *  W  +  arc  tan 


/  3  -  m  'N 

V  t/pe  -  s  +  |y|  / 


(39) 


4ama 

E  Vs2  -  y2 


(40) 


m 


( 


<1o4 


8a 

«Mq 


+  y  - 
-  y 


s 


+  rnMoP© 


-  y  - 

s  +  y 


s 


+ 


arc  tan 


/ 


s  -  y 


t/3( 
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+  y 


/  3  +  y  _  21  ^ 

J  t/pft  -  y  -  s  2  J 


(42) 


The  load  distribution  across  any  section  is  given  by  equations  (39) y 
(4l),  and  (42)  from  the  time  t/pe  =  0  to  (t/pe)i,  where  the  term 
(t/p0)1  is  equal  to  2s  or  s/m(Mo+l)p0,  whichever  is  smaller. 
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(At  t/(Be  =  2s  the  span  wise  section  has  just  reached  the  secondary 
waves  shown  in  figure  6,  and  at  t/(3e  =  s/mpe  (Mq+1  )  the  spanwise  section 
has  just  reached  the  spherical  wave  which  started  from  the  apex. )  From 
(t/fBg^  to  (t/pe)2  =  s/mPe(M0- 1),  the  loading  has  not  been  determined, 

and  from  t/pe  =  (t/pe)2  to  t=  co 
the  loading  is  the  steady— state 
value  given  by  equation  (40). 

The  sketch  shows  this  initial  and 
final  load  variation  plotted  as  a 
1  function  of  the  parameter  t/pe. 

At  the  beginning  of  the  motion 
the  loading  is  constant  across 
the  span,  but  this  type  of  dis¬ 
tribution  is  quickly  modified  and 
the  shape  of  the  curve  tends 
toward  the  "inverted  elliptic" 
loading  given  by  equation  (40) 
and  shown  in  the  sketch  as  the 
distribution  at  t/pe  =  (t/pe)2. 

In  fact,  when  the  span  has 
traveled  a  distance  such  that 
t/pe=2s,  the  expression  for  the 
/  loading  given  by  equation  (42) 
becomes 

(%')  16  {h3) 
Vlo/JL  =  2s  it  s2— y s 

Pe 

which  differs  from  the  value  given  by  equation  (4-0)  only  by  a  constant 
of  proportionality.  Both  before  and  after  the  time  t/pe=2s  the  shape 
of  the  loading  curve  varies  from  the  simple  type  represented  by  equa¬ 
tion  (43),  but  the  trend,  and  to  a  certain  extent  the  rapidity  of  the 
trend,  is  clearly  established* 

The  average  spanwise  loading  P0,  introduced  by  equation  (20a),  can 
now  be  determined  for  certain  regions*  Hence,  if  the  notation 


=  t/s 


2sa  sLg 


is  adopted,  there  results  for  the  early  part  of  the  motion,  that  is,  for 

0  <  r/p  <  (  r/Pgh 


=  2 


^  J  +  4mpeMQ  (  T/pe) 
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Equation  (45)  was  derived  by  integrating  equations  (39),  (4l),  and  (42). 


For  values  of  t  /pe  >  l/pem(M0-l ) 
Hence 


equation  (40)  is  valid. 


P 


o 


2JtmMp 

E 


(46) 


The  sketch  indicates  the 
magnitude  of  this  average  load 
for  loth  large  and  small  values 
of  T/pe .  Notice  that  for  small 
values  of  t/Pq  it  is  suffi¬ 
cient  for  the  establishment  of 
the  curve  to  specify  the  parameter 
mMo0e,  but  for  large  values  an 
additional  parameter  must  be  given 
(such  as  Mq  in  the  sketch). 
Notice,  further,  that  in  spite  of 
the  large  variation  in  the  distri¬ 
bution  of  the  loading,  as  shown  in 
the  previous  sketch,  the  average 
value  PQ  varies  linearly  through¬ 
out  the  intervals  considered.  This 
result  is  similar  to  the  one 
obtained  for  triangular  wings  with 
supersonic  edges  and  is  given  in 
equations  (24). 


INDICIA!  LOADING  ON  VERY  SLENDER  TRIANGULAR  WINGS 


Analysis 


In  the  first  section  of  this  report  entitled  "The  Problem  and  the 
Nature  of  its  Solution,"  it  was  pointed  out  that  if  the  wing  is  slender 
(i.e.,  has  a  small  ratio  of  span  to  chordwise  length)  the  ba3ic  partial 
differential  equation  (l)  can  be  approximated  by  the  equation  which  was 
previously  introduced  as  equation  (4),  thus: 


9tt  ^yy  'Pzz  ® 


The  boundary  conditions  appropriate  to  this  problem  will  now  be  examined 
in  some  detail. 
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Just  as  in  the  previous  sections  of  this  report,  consider  a  triangu¬ 
lar  wing  which  is  at  rest  for  t  <  0,  starts  suddenly  to  move  at  a  for¬ 
ward  velocity  equal  to  VQ  at  t=0,  and  continues  at  this  same  velocity 
for  t  >  0.  It  should  he  emphasized  that  in  this  case,  VG  te 

subsonic  or  supersonic.  A  section  in  the  spanwise  direction,  as  for 

instance  section  AA  in  the  sketch, 
projects  into  the  yt  plane  as  a 
>1\.  „„„  tXn  m  narrow  rectangular  strip  along  the  t 

/  V  °r  axis.  Since  equation  (4)  has  been 

^ _  _  /  7_\_  _  _ A  derived  on  the  assumption  that  the 

71  \  velocity  gradients  in  the  y,  z,  and 

/  |  |  \  t  directions  are  independent  of  the 

/  \  gradient  in  the  x  direction,  the 

/  \  \  A  boundary  conditions  along  the  strip 

|  I  shown  in  the  sketch  are  independent 

s=rr)X  of  those  on  other  strips  projected. 

I  'fjK  I  from  spanwise  sections  along  the  wing. 

I  I  Hence,  the  problem  is  to  find  a  solu— 

) I  ^  ^ _  y  tion  to  equation  (4)  which  will  make 

■j-r  y  i  1  ^  q>  constant  over  the  strip  and  at  the 

'//'//  same  time  will  satisfy  the  other  con- 

yv  f //  ditions  listed  under  equation  (l). 

'/////  In  the  lifting-surface  analog  this. 

/ / /  / / /  corresponds  to  the  problem  of  finding 

// >  / //  the  velocity  potential  over  a  flat 

/////j  rectangular  wing  of  low  aspect  ratio 

/////  situated  in  a  free  stream  moving  at  a 

'  Mach  number  equal  to  «/"5*  Solutions 

1 '  to  the  latter  problem  can  be  obtained 

by  various  techniques,  and  so  the  pro¬ 
cedure  will  be  first,  to  find  the  potential  for  the  steady-^tate,  flat, 
rectangular  wing,  and  then,  by  analogy,  to  convert  this  to  the  solution 
for  the  very  slender  triangular  wing  in  unsteady  motion. 

The  steady-state,  lifting-surface  problem.-  Lifting-surface  solu- 
tions  for  the  loading  on  a  rectangular  wing  traveling  at  supersonic 
speeds  have  been  developed  for  regions  1,  2,  and  3  of  ^ 

Busemann  and  others),  and  by  means  of  these  solutions  the  load  dis  r 
button  on  a  spanwise  section  of  the  triangular  wing  can  e  e  ermine 
to  a  time  necessary  for  sound  to  travel  that  span  length.  s, 

however,  the  solution  becomes  considerably  complicated  by  the  increasing 
number  of  reflections  from  the  edges.  Keference  8. gives  solutions  for 
the  loading  on  a  rectangular  wing  in  region  4  and  indicates  methods. for 
extending  the  solution  to  regions  farther  along  the  wing.  Alrea  y  m 
region  4,  however,  the  expression  is  cumbersome  and  in  higher-numbered 
regions  the  expressions  become  difficult  to  manipulate.  These  methods, 
therefore,  will  be  discarded  in  favor  of  a  more  approximate  but  simpler 
analysis. 


If 

span^  and 


is  the  distance  along  the  chord,  y  the  distance  along  the 
s  the  semispan,  then  the  solution  for  regions  1,  2,  and  3  °f 
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figure  8  can  te  written  (for  convenience,  see  reference  5): 

Region  1 

^  =  4a  (47a 

^-o 

Region  2 

=  §1  arc  tan  /ZZZEE  (47b 

<Lo  *  v  x  —  s  +  |  y  | 

Region  3 

—  =  —  (arc  tan  [  8+y  +  arc  tan  /  ^  (Vfc 

qQ  n  \  V  x-^3-y  v  x-s+y  2  / 

As  x  increases  (i.e.,  for  higher— numbered  regions  in  figure  8) 
it  is  reasonable  to  assume  that  the  spanwise  variation  of  loading  is 
relatively  unimportant  —  except  that  it  be  Tf smooth"  and  fall  to  zero  at 
the  side  edges  —  and  the  chordvise  variation  of  loading  is  dominant. 
Assume,  therefore,  that  the  loading  is  given  by  the  relation 
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Spanvise,  this  has  the  variation  shown  in  the  preceding  sketch;  chord- 
wise  it  is  as  yet  arbitrary.  To  fix  the  chordwise  distribution  the 
value  of  f(x/s)  will  be  determined  so  that  the  vertical  induced  velo¬ 
city  along  the  center  line  is  constant  and  equal  to  — ' V0a. 


The  solution  to  this  somewhat  artificial  problem  approaches  the 
exact  solution  to  the  steady-state  lifting-surface  problem  for  a  flat 
rectangular  wing  along  sections  far  behind  the  leading  edge;  closer  to 
the  leading  edge  it  only  approximates  the  exact  solution;  and,  of  course, 
in  the  vicinity  of  the  leading  edge  it  will  be  least  representative. 

But,  on  the  other  hand,  the  exact  solution  is  known  in  the  vicinity  of 
the  leading  edge  and  it  turns  out  that  the  solution  of  the  problem  posed 
above  forms  a  reasonable  continuation  over  the  remainder  of  the  wing. 

The  velocity  potential  for  the  problem  which  has  been  set  can  be 
readily  expressed  in  terms  of  an  integration  of  elementary  horseshoe 
vortices  over  the  plan  form.  Since  the  Mach  number  equals  J 2,  then 
according  to  reference  9 > 

(x-xi)  (Ap/qQ)  dxidyi 
(yi2+z2)  J (x-xi)2  -y 12— z2 


where  A  is  the  area  on  the  wing  within  the  forecone  from  the  point 
P (x,y, z ) ,  at  which  cp  is  to  be  determined  (the  shaded  area  in  the 
sketch). 

The  simplification  of  the  last  expression  is  given  in  reference  9* 
The  result  is  the  integral  equation 

1  =  f (q )  +  |  tin  )  Gr  (r1-n1)dri1  (49) 

J0 

where  t)=x/s  and  G  is  given  by 

1  < 
ft* 

q-ni  <  1 

i 

ki  =  5  ks  = 


O(Mi)  = 


Ei 

E2-(l-k2; 

k2 


The  modulus  of  E1  is  ki  and  the  modulus  of  Kg  and  E2  is  k2. 
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The  solution  of  equation  (49)  for  f (-q )  is  not  difficult  when 
numerical  methods  are  used.  For  intervals  of  tj ^  equal  to  0.2,  the 
result  is  given  in  tabular  form  in  appendix  A,  and  also  in  the  sketch. 


This  result  can  be  improved  in  the  interval  0  <  q  <  2  by  means  of  the 
correct  solution  in  that  interval  given  by  equations  (47).  Comparison 
of  the  spanvise  average  of  .the  loading  given  by  equation  (48)  with  that 
which  can  be  derived  from  equations  (4-7)  gives  an  equivalent  f(q)  in 
the  interval  which,  when  used  in  equation  (48),  will  give  the  correct 
value  of  the  average  span  loading.  The  sketch  also  shows  a  curve  for 
this  equivalent  f(q)  which  starts  at  k/n  and  falls  linearly  to  zero 
at  q=2. 

By  using  the  sketch,  or  the  results  listed  in  appendix  A,  the 
loading  over  a  low-aspect— ratio  rectangular  wing  flying  at  a  Mach  number 
equal  to  can  be  estimated.  Of  particular  interest  is  the  damped 

oscillatory  nature  of  the  load,  falling  to  zero  at  one  span  length 
behind  the  leading  edge  and  taking  alternately  negative  and  positive 
values  beyond  this  point.5  A  somewhat  different  approach  to  this 
problem  (reference  11)  has  recently  led  to  a  solution  very  like  the  one 
given  here. 

The  unsteady  analog,  sinking  wing.—  The  first  step  in  deriving  the 
unsteady-flow  results  for  the  sinking  wing  from  the  steady  solution  is 
to  replace  x  with  t.  In  equation  (49)  this  corresponds  to  replacing 
q  with  t  where  t  is  equal  to  t/s  (equation  (44).  The  second  step 
is  to  rederive  the  expression  for  loading  coefficient  since  in  the 
time-varying  problem  it  is  expressed  in  a  somewhat  different  manner  than 
in  the  steady— state  analog.  In  the  unsteady  case,  as  the  triangular 
wing  moves  through  a  fixed  reference  plane  the  local  span  intersecting 


A  discussion  of  this  aspect  of  these  results  is  contained  in  refer¬ 
ence  10. 
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this  plane  grows  as  a  function  of  time  and.  equation  (2)*  which  repre¬ 
sents  the  partial  derivative  with  respect  to  time  with  x  fixed*  must 
he  expanded  to  the  form 


Ap  =  2 

4o  VqMo 


8Acp 

2 

/ 

[ 

'8a$ 

+ 

8a<p 

|_at  J 

x  voMo 

V 

_  ^ 

S 

8s 

where 


and 

mM0* 


8ACp 

and 

8Aq? 

_  8t  _ 

s 

8s  _ 

t*  respectively.  Since 
and  there  results 


indicate  derivatives  taken  at  constant  s 
3  is  equal  to  m(x  +Mot)*  8s /8t  equals 


Ap  2 
4o  '  V°Mo 


( 


rsAcp- 

[8t  j 


+  mMo 


SAcp 

ds 


) 


(50) 


8a<P~ 
8s_ . 


In  the  steady-state  problems  an  analog  to  the  term  involving 

is  missing*  and  the  loading  coefficient  is  given  entirely  by  an  opera- 

8A<P 

-  s 


tion  equivalent  to 


8t 


It  is  necessary*  therefore*  to  operate 


further  on  the  solution  given  for  the  loading  in  the  steady— state 
problem  to  obtain  the  solution  for  the  loading  in  the  unsteady  problem. 
But 


8A(p 

,A.r\ 

8Acp 

_8s  _ 

t 

8s  J0 

.  8t  _ 

dti 


so  that  if  the  notation 


is  adopted  (where 


,1o  / 


represents  the  loading  in  the  analogous 


steady-state  problem)*  then  the  expression  for  the  unsteady  loading  can 

/.  \ 

be  gi  ,ren  in  terms  of  (  —  )  by  the  equation 

vio/a 


Ap  _  1 

q.o  Mo 


V\ 

V. 


+  mMo 


ds 


t  '&)  dt 

va 


(51) 
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By  the  application  of  equation  (51)  to  equations  (47  )j  the  loading 
for  the  various  regions  of  the  unsteady  wing  can  he  found.  For  region  2 
there  results 


Ap 


_1 

M 


o  L 


8a 


arc  tan 


n 


J 


t— s  +  J  y ;  os 


s— |y| 


4a  dti  + 


arc 

J t 


tan  J — — 421. 


ti-s+  I  y  I 


dt 


which  becomes 


Ap  _  8a 

q0  “  «M0 


f  t — S- 

rv^r 


t — S+  1  y  | 


1  yi 


+  arc  tan 


J- 


S-  I  y  I 


t-a+  |y [ 


The  loading  coefficient  can  he  similarly  derived  in  the  other  regions  so 
that  finally 


Region  1 

Ap  4a 

lo  Mo 


Region  2 


Ap  _  8a 

~*“Mo 


t~^+  |y 


+  arc  tan 


42= - 

V  t— s+ 


ki_ 

|y| 


) 


Region  3 


4p 

<lo 


(  mMo  /EEZ  +  mM0  /EE*  ♦ 

*M0\^  s-y  o  ^  s+y 

arc  tan  +  arc  tan  [  —  —  ''j 

V  t-s+y  V  t— s— y  2  / 


(52a) 


(52b) 


(52c) 
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For  the  interval  t  >2s  equation  (48)  must  he  considered.  By  means 
of  equation  (51),  the  expression  for  the  loading  coefficient  can  he 
written 


A. (z)  , «<, I- /"V.  f  &,')  A. (z) 

lo  Mo  L  V  W  dsj  Vs  /v  W 


dt 


which  becomes 


Ap  =  4a 
do  Mo 


(l-mMo  t) 


'1- 


f(r)  + 


mMo 


0 


y1  f(Ti)dT1l  (53) 


where  f("r)  is  the  solution  to  equation  (49).  Notice  that  for  large  T 
(when  the  loading  has  reached  its  steady  state), f(f)  is  zero  and 


L 


f(Ti)dT1  is  unity, 
hy  the  equation 


(See  appendix  A.)  Hence,  the  loading  is  given 


Ap  _  4oms 
<10  J  S2_y2 


which  is  the  steady— state  value  for  a  slender  triangular  wing  (equa¬ 
tion  (40)  when  E  =  l). 

It  is  now  possible  to  derive  the  average  span  loading  PD  as 
defined  hy  equation  (20a),  thus 


Po 


Placing  equations  (52)  and  (53)  in  this  expression,  it  is  found  that 
for  0  <t<  2 

Po  =  2(2— t)  +  4mMoT  (54a) 


and  for  v  >  2 


PD  =  Jt(l-nM0T)f(r)+2  mMoit  P  f(Tx)  dTX  (54b) 

J0 

Since  the  values  of  PQ  given  hy  equations  (54a)  and  (54b)  were 
derived  using  different  methods,  their  magnitudes  at  t  =  2  are  not 
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equal.  The  final  curve  for  Po  must  be  constructed  by  fairing  the  solu¬ 
tion  for  t^2  into  that  for  t  >2.  The  accompanying  sketch  shows 
these  re suit 3  together  with  the  final  curve  chosen  (solid  line). 


The  unsteady  analog,  pitching  wing.—  When  the  wing  is  pitching  at  a 
steady  rate  about  its  apex,  the  equation  for  the  vertical  induced  veloc¬ 
ity  on  the  plan  form  is 

wu  =  — (x  +  Mot)0 


so  that  the  a  in  the  steady-state  equations  (47)  and  (48)  must  be 
replaced  by  0s/mVo.  Since  the  loading  coefficient  is  still  given  by 
equation  (51),  there  results  for  the  conversion  of  equation  (48)  the 
expression 


*0 


49 


mM„Y 


0  0  L- 


</ s2— y2  f ( t)  +  mMo 


a_ 

8s 


f(Ti) 


dTi 


and  this  can  be  reduced  to  the  form 
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(1-mMoT)  Zs2-y2  f (t)  + 

q_  mMoVo 


m  z'28^  r  f(T1)  dT1  (55) 

v°  \V  a  2— y2/ 

As  in  the  discussion  of  equation  (53),  it  can  be  seen  that  equation  (55) 
becomes  for  the  steady  state  (T  large) 

Ap  _  40  ^ 2sg-y 2  ^ 

^  =  ^o 

and  this  can  be  shown  to  agree  with  the  steady-state  slender-wing  results 
given  in  reference  12 • 

It  is  now  possible  to  derive  the  average  span  loading  Pi  as 
defined  by  equation  (20b) 


pi  =  ;so  r^iy 

2s20  J— s  ^o 

Using  equation  (55),  one  finds  for  t  >2 


% 


it(l-mM0T) 


f  (  t)  +  3mMoit 


dTi 


(56a) 


and  a  similar  analysis  based  on  equation  (47)  yields 
for  0<  t  <  2 

I3!  =  —■  (  2— T+4mMQT—  i  mMoT2)  (56b) 

As  in  the  case  for  PQ,  the  two  equations  for  Px  do  not  join  at  t=2 
and  the  final  curve  must  be  constructed  by  fairing  the  solution  for  t<,2 
into  that  for  t  >  2. 
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Discussion  of  Results 

It  is  now  possible  to 
assess  the  accuracy  of  the 
solution  for  very  slender 
wings  in  the  interval  0  <;t  £2 
hy  comparing  equations  (52) 
and  (54a)  with  equations  (hi), 

(39),  (42),  and  (45),  the 
exact  solutions  for  this  region 
derived  in  the  preceding  sec¬ 
tion.  It  is  apparent  that  the 
approximate  solution  differs 
from  the  exact  only  by  a 
stretching  factor  in  the  t 
direction.  Hence,  if  t  is 
replaced  by  r/pe  and  m 
(note  m  is  proportional  to 
y/t)  by  mj3e,  where  Pe  is 
given  by  equation  (38),  then 
equations  (52a),  (b),  and  (c) 
are  identical  with  equa¬ 
tions  (4l),  (39),  and  (42),  respectively,  and,  of  course,  equation  (54a) 
corresponds  to  equation  (45). 

This  rather  remarkable  result  can  be  enlarged  upon  from  another 
viewpoint.  Suppose  that  in  the  steady-state  analog  problem  the  wing  had 
been  flying  at  some  Mach  number  other  than  *T 2,  say  M©.  The  solution 
to  such  a  new  problem  could  be  obtained  from  the  old  one  merely  by 
applying  the  Prandtl-Glauert  correction,  that  is,  by  stretching  all  dis¬ 
tances  in  the  x  direction  by  the  factor  l/pe  where  Pe2=  |l“M©2|. 
Such  a  procedure  would  convert,  for  example,  equation  (54a)  to  the  form 

Po0e  =  2(2—  t  /pQ )  +  4m0eMo  T/0e 

Finally,  if  P0  ia  adjusted  so  that  P0=4  at  t  =0,  there  results 

PQ  =  2(2-  r/0e)  +  4mMoT 

which  is  exactly  the  answer  given  by  equation  (45).  It  is  possible  to 
simplify  the  statement  of  this  procedure  by  simply  remarking:  The  exact 
results  for  Ap/q0  or  PQ  in  the  interval  0  <  T/pe  <  2  can  be 
obtained  from  the  approximate  results  for  a  very  slender  wing  by  making 
an  effective  Mach  number  correction  to  the  right-hand  side  of  equa¬ 
tions  (52)  or  (54a),  respectively. 
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It  is  interesting  to  pursue  this  concept  even  further.  Consider  a 
spanwise  section  of  a  triangular  wing  as  time  increases  from  the  starting 

impulse.  The  primary- 
wave  fronts  emanating 
from  either  side  pass 

Actual  trace  of  acroag  the  section 

primary  wave  front  forming  the  Mach  lines 

in  the  steady— state 

L _  j  _ ,  rectangular— wing 

_ _  %  analogy.  For  very 

\  /  '  y/s  slender  wings  these 

\  /V  lines  make  a  45°  angle 

VS.  J  with  the  trace  of  the 

side  edge  and  are  used 

\  /,  to  divide  the  plan 

'  form  into  regions  as 

/  *  in  the  sketch.  Now 

y'  ' '  3“  find  the  actual  posi— 

/  /  \/\ye  tion  of  these  primary 

/  y  /  n.  wave  fronts  as  they 

/  /  ^  \  form  a  trace  on  the 

'  /  |  \  \  I  section  in  the  yT 

/  \  plane.  A  straight- 

Trace  Of  primary  \  forward  calculation 

wave  front  obtained  x  shows  that  these  lines 

from  s/ender-wmg  j  equal  to  arc  tan  i/pe 

theory  1  with  the  trace  of  the 

I  side  edges.  Hence  the 

jy.  effective  Mach  number 

which  is  used  to 
correct  the  slende ra¬ 
wing  results  in  the 

interval  0  <  t  <  2f3e  is  that  which  makes  the  Mach  lines  of  the  steady— 
state  analogy  coincide  with  the  actual  trace  of  the  primary  wave  fronts. 


//  45° 

tori’! L 
\  fie 


Trace  of  primary  \ 
wave  front  obtained  \ 
from  slender-wing 
theory  \ 


INDICIAL  LIFT  AND  PITCHING  MOMENT  ON 
VERY  SLENDER  TRIANGULAR  WINGS 

Analysis 


The  lift  coefficient  for  the  sinking  wing  is  given  in  the  notation 
introduced  in  equation  (21a)  "by  the  equation 


Pod* 
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where  PQ  has  been  determined  in  the  last  section  as  a  function  of 
t =t /s »  Consider  the  situation  at  a  certain  fixed  time  and  let  the  x 
coordinate  in  the  above  formula  be  fixed  in  the  wing.  Then  set 


and  as  before 


(57) 


where  tq  is  the  number  of  wing-chord  lengths  traveled.  In  this  way 
the  equation  for  lift  coefficient  becomes 


and  by  a  similar  analysis  the  pitching-moment  coefficient  taken  about 
the  apex  can  be  written 


(60) 


The  equation  for  the  lift  and  pitching-moinent  coefficients  (where 
again  the  pitching  moment  is  taken  about  the  apex)  on  a  pitching  wing 
are 


and 


(62) 


Discussion  of  Results 


The  values  of  PQ  and  Pi  were  taken  from  curves  similar  to  the 
sketches  in  the  last  section  (using  the  faired  curves  in  the  vicinity  of 
t=2)  and  the  results  for  the  indicial  lift  and  pitching  monent  in  terms 
of  tq,  the  number  of  chord  lengths  traveled,  are  shown  in  figures  (9) 
and  (10)  for  a  value  of  mMo  equal  to  l/8.  The  results  are  all 
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APPENDIX  A 

SOLUTION  OF  AN  INTEGRAL  EQUATION 


For  convenience  in  applying  equations  (54b)  and  (56a)  of  the  text, 
a  table  of  values  of  the  function  f ( r) ) ,  obtained  by  numerical  solution 
of  equation  (4-9),  is  given  here.  The  values  of  the  integral, 

fQ  f(Hi)cLTji,  are  also  listed. 


I 

f(n) 

/0Vil  i)ani 

n 

f(n) 

J^fCniHni 

0.0 

1.0000 

0.0000 

4.0 

-0.1307 

1.0124 

.2 

.9899 

.1990 

4.2 

-.1008 

.9893 

.4 

.9597 

.3940 

4.4 

-.0716 

.9720 

.6 

.9087 

.5808 

4.6 

-.0446 

.9604 

.8 

.8356 

•  7552 

4.8 

-.0207 

.9539 

1.0 

.7339 

.9122 

5.0 

-.0008 

.9517 

IB 

.6032 

1.0459 

5-2 

.0149 

.9531 

1.4 

.4597 

1.1522 

5.4 

.0265 

.9573 

HEU 

.3188' 

1.2300 

5.6 

.0340 

.9633 

1.8 

.1880 

1.2807 

5.8 

.0379 

.9705 

2.0 

.0724 

1.3067 

6.0 

.0388 

.9782 

2.2 

-.0245 

1.3115 

6.2 

.0373 

.9858 

2.4 

-.1008 

1.2990 

6.4 

.0338 

.9929 

2.6 

-.1562 

1.2733 

6.6 

’.0292 

.9992 

2.8 

-.1919 

1.2385 

6.8 

.0237 

1.0045 

3.0 

-.2098 

1.1983 

7.0 

.0180 

I.OO87 

1 

-.2126 

1.1561 

7-2 

.0124 

I.0117 

3.4 

-.2031 

1.1145 

7.4 

.0072 

1.0137 

3.6 

-.1843 

1.0758 

7.6 

.0026 

1.0146 

3.8 

-•1523, 

i.04l4 

7.8 

-.0012 

1.0148 
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wing  with  supersonic  leading  edges. 
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Figure  4.-  indicia /  aerodynamic  characteristics  of 
sinking  wings  with  supersonic  leading  edges.  M0-  2. 
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(b)  Pitching  moment  about 
leading  edge  or  apex . 


Figure  4.  -  Continued . 
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(c )  Location  of  center  of  pressure 
with  reference  to  the  leading 
edge  or  apex. 


Figure  4.-  Concluded. 


Chord  lengths  traveled,  tc 
(a)  Lift. 

Figure  5.-  indicia l  aerodynamic  characteristics  of 

wings  with  supersonic  edges  pitching  about  lead¬ 
ing  edge  or  apex.  M«s  2. 
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(b)  Pitching  moment  about 
lead ing  edge  or  apex. 


Figure  5.-  Continued. 


■4  .8  L2  1.6 

Chord  lengths  traveled ,  r0 

(c)  Location  of  center  of  pressure 
with  reference  to  the  leading 
edge  or  apex. 


5.-  Concluded. 
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triangular  wing  with  subsonic  leading  edges. 


y-(x+  M0t)//5 
y=m(x+  M0t) 


Regions  used  in  the  discuss i> 
-ratio  rectangular  wing. 
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Chord  lengths  traveled,  rB 


fa)  Lift. 


(b)  Pitching  moment  about  apex. 


Figure  9—  indicia i  aerodynamic  characteristics  of  sinking 
triangular  wings  with  slender  plan  forms .  mM>-  i/8. 
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Figure  iO.~  Indicia!  aerodynamic  characteristics  of  tri¬ 
angular  wings  with  slender  plan  forms  pitching  about 

apex,  m Mo  -  /  /8. 
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